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Abstract. - We investigate the distribution of the resonance widths V(T) and Wigner delay 
times V(tw) for scattering from two-dimensional systems in the diffusive regime. We obtain the 
forms of these distributions (log-normal for large tw and small V, and power law in the opposite 
case) for different symmetry classes and show that they are determined by the underlying 
diffusive classical dynamics. Our theoretical arguments are supported by extensive numerical 
calculations. 



Quantum scattering has been a subject of intensive research activity both in mesoscopic 
physics and in Quantum Chaos in the last years [1-7] . Among the most interesting quantities 
for the description of a scattering process are the Wigner delay times and resonance widths. 
The former quantity captures the time-dependent aspects of quantum scattering. It can 
be interpreted as the typical time an almost monochromatic wave packet remains in the 
interaction region. It is related to the energy derivative of the total phase shift &{E) = 
— i lndet S(E) of the scattering matrix S(E), i.e. tw(E) — ■ Resonances are defined as 

poles of the ^-matrix occurring at complex energies £ n — E n — |r n , where E n is the position 
and r n the width of the resonance. They correspond to "eigenstates" of the open system that 
decay in time due to the coupling to the " outside world" . 

For chaotic/ballistic systems Random Matrix Theory (RMT) is applicable, and the dis- 
tributions of resonance widths V(T) and Wigner delay times Vijw) are known [2]. As the 
disorder increases the system becomes diffusive and the deviations from RMT become increas- 
ingly apparent. In the strongly disordered limit where localization dominates, the distribution 
of resonances ViT) [4] and delay times V(jw) [5] were found recently. At the same time, an 
attempt to understand systems at critical conditions, was undertaken in [6,7]. For diffusive 
mesoscopic samples, however, there is no study of V(T) and V(tw) besides Ref. [3] where 
the authors have focused on the tails of V(T) for a quasi-lD system in the diffusive regime. 
This study is important for diffusive random lasers, where the knowledge of short resonance 
width distribution determines the properties of lasing thresholds [8], as well as for various 
other applications like mesoscopic capacitors [9], microwave cavities [10] and chaotic optical 
cavities [11] where most of the theoretical treatment is limited by RMT. We point here that 
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current developments of microwave experiments in random dielectric media in the diffusive 
regime [12] may allow us a direct comparison between theory and experiment. Finally we note 
that an analogous theoretical study for the distribution of wavefunction intensities of closed 
systems was put forward very early in the development of the mesoscopic physics [13] and 
recently their statistical properties have become clear [14-18]. 

In this paper, we study V(T) and V(tw) for 2D open systems in the diffusive regime. We 
will show that these distributions are determined by the diffusive classical dynamics of the 
corresponding closed system and depend on the time-reversal symmetry (TRS). Specifically, 
the resonance width distribution V(T) is given by 

V(T<T c i) ~ exp(-C /3 (lnr) 2 ) , where C fj ~f3D 

^(r>r ci) ~ (!) 

where D is the classical diffusion coefficient (which is proportional to the Thouless conductance 
for disordered 2D systems), L is the linear size of the system, r c ; is the classical decay rate [19] 
which is inversely proportional to the Thouless time L 2 /D and /3 denotes the symmetry class 
[j3 = 1(2) for preserved (broken) TRS]. For the distribution of Wigner delay times we obtain 

'PiTW <T a V ) ryj exp(-0-/T W ) 

T(T W >r- l 1 ) ~ expt-C^lnny) 2 ) (2) 

where a is some constant of order unity. Our results (11121) are applicable both for chaotic and 
disordered systems, provided that: (i) transport between scattering events may be treated 
semiclassically, (ii) a particle scatters many times while it traverses the system and (iii) strong 
localization effects can be neglected. 

To numerically demonstrate the above predictions we use a model of the 2D Kicked Rotor 
(KR) with absorbing boundary conditions. The corresponding closed system is described by 
the time-dependent Hamiltonian [17,20] 

2 

H = H + kV^28(t-m) ,if ({A}) = Ef £ i ( 3 ) 

m i—1 

V{{9i}) = cos(0i) cos(6> 2 ) cos(a) + - sin(26>i) cos(26» 2 ) sin(a) 

where Ci denotes the angular momentum and 0j the conjugate angle of one rotor. Here the 
kick period is one, k is the kicking strength, while pi is a constant inversely proportional to 
the moment of inertia of the rotor. The parameter a breaks TRS. The motion generated by 
(|3*|) is classically chaotic and for a sufficiently strong kicking strength k there is diffusion in 
momentum space with diffusion coefficient D = lim t ->oo < £ 2 (t) > ft — \ [20]. 

Quantum mechanically, this system can be described by a finite-dimensional evolution 
operator for one period 

U = e (-iHo({£*})/2) e (-iV({<>*})) e (-iHo({£i})/2) ( 4 ) 

where we put % = 1. The eigenvalues of U are A„ = e lu)n and lie on the unit circle; u n are 
known as quasi-energies and are dimensionless. The corresponding mean quasi-energy spacing 
is A = 2ir/L 2 , where L is the linear size of the system. A detailed presentation of our model 
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can be found in [17], where the statistical properties of wavefunction intensities in the regime 
of interest were investigated and found to be in agreement with the theoretical predictions for 
mesoscopic diffusive samples [14,15]. Therefore the 2D KR is a representative model both for 
2D disordered and chaotic systems. 

In this paper, we turn the closed 2D KR model into an open one. To this end we impose 
absorption at the boundary of a square sample of size L x L in the momentum space. In other 
words, every time that one of the components of the two dimensional momentum (£1,^2) 
takes on the value 1 or L, the particle is absorbed without coming back to the sample. Using 
a recently proposed recipe [21] we can write down the corresponding scattering matrix S in 
the form [22] 

S(u) = -WUe iu - 1 —— W f ,P^I-W^W (5) 

where / is the L 2 x L 2 unit matrix and W is a M x L 2 matrix. It has only M non-zero 
elements which are equal to one and describe at which "site" of the L x L sample we attach 
M "leads" [in our case M — 4(L — 1)]. Here W^W is a projection operator onto the boundary, 
while P is the complementary projection operator. The scattering matrix Sij given by Eq. lf5"jl 
can be interpreted in the following way: once a wave enters the sample, it undergoes multiple 
scattering induced by [/ — e lw -PC/] -1 = J2^=o (e luJ PU) until it is transmitted out. It is 
clear therefore that the matrix U — PU propagates the wave inside the sample. However, 
contrary to the closed system in which the evolution operator is unitary, the absorption breaks 
the unitarity of the evolution matrix U so that all eigenvalues A move inside the unit circle. 
Therefore each eigenvalue can be written in the form A„ = e tuJn = exp(—ioj n — T n /2) where 
r„ > is the dimensionless resonance width of an eigenstate. 

The Wigner delay time can be expressed as the sum of proper delay times r q . The latter 
are the eigenvalues of the Wigner-Smith operator written in our case as [22] 

1 d S* 

Q(u>) = -S^— = -e~ iu WK^W^WKU^KW^ (6) 
1 dui 

where K = (P - We'™) ~\ 

Below we present our theoretical considerations and compare them with the numerical 
data obtained for the 2D KR model. The parameters of the model were chosen in such a 
way that the conditions (i)-(iii) discussed above were fulfilled (see [17]). In order to improve 
our statistics, we randomized the phases of the kinetic term of the evolution operator 0} and 
used a number of different realizations. In all cases we had at least 60000 data for statistical 
processing. 

Resonance widths distribution. We start our analysis with the study of resonance width 
distribution V(T) for T < T c j. The small resonances V < A can be associated, with the 
existence of pre- localized states of the "isolated" system. The latter states show strong lo- 
calization features in contrast to the typical states in the diffusive regime and are considered 
as precursors of localization. They consist of a short-scale bump (where most of the norm is 
concentrated) and they decay rapidly in a power law from their center of localization [14, 15]. 
One then expects that states of this type with localization centers at the bulk of the sam- 
ple are affected very weakly by the opening of the system at the boundaries. In first order 
perturbation theory, considering the opening as a small perturbation we obtain 
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FIG.l: (a) The distribution of resonance widths (plotted as V(l/T) vs. 1/r) for T < r c ; for two representative 
values of D. The system size in all cases is L = 80. Filled symbols correspond to broken TRS. The solid lines 
are the best fit of Eq. for /3 = 1(2) to the numerical data, (b) Coefficients On vs. D. The solid lines are 
the best fits to Gp = ApD + Bp for = 1(2). The ratio R = A 2 /A 1 = 1.95 ± 0.03 



where |^]/(L)| 2 is the wavefunction intensity of a pre-localized state at the boundary. At the 
same time the distribution of 9 = l/VLf(L) for large values of the argument is found to be of 
log-normal type i.e. V{9) ~ exp (—ir 2 Dln 2 (9 2 )) [15]. Using this together with Eq. (EJ we get 
V(l/T) - exp (-ir 2 Din 2 We would like to stress that the expression for V(9), must 
be corrected by including the TRS factor (3 in the exponent. This is due to the fact that the 
Optimal Fluctuation Method, which was used to derive the above expression for P(9), does 
not describe the effect of breaking TRS in a correct way [17, 18]. Taking all the above into 
account we end up with the expression given in Eq. p|l. 

The numerical data reported in Fig. 1 support the validity of the above considerations. 
However, we mention that the perturbative argument is valid only for the case of very small 
resonances i.e. T < A, whereas our numerical data indicate that one can extend the log-normal 
behaviour of V(T) up to resonances with A < T < T c i. 
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FIG. 2: (a) Resonance width distribution V(T) for preserved TRS and D = 20.3 (o) and D = 33.5 (0)- The 
corresponding filled symbols represent V(T) for broken TRS and the same values of D. The dashed (solid) 
vertical line mark the classical decay rate r c ; [19] for D = 20.3(33.5). (b) The "Pi n t(r) for a sample with 
nine leads (lower curve). For comparison we plot also the "Pi nt (r) for the same sample but when we open the 
system from the boundaries. The dashed lines correspond to the theoretical predictions J5J and 1101 . 



Next we turn to the analysis of V{Y) for Y > Y c i. In this case, the resonances are overlap- 
ping, the corresponding eigenstates are strongly non-orthogonal and first-order perturbation 
theory breaks down [23]. In Fig. 2a we report our numerical results for V(Y) with preserved 
(broken) TRS for two representative values of D. An inverse power law V(Y) ~ r~ 15 is evi- 
dent in accordance with Eq. (JTJ [24] . From Fig. 2a it is clear that this part of the distribution 
is independent of the symmetry class, in contrast to the small resonance distribution. 

The following argument provides some understanding of the behaviour of V(Y) for Y > r c ;. 
First we need to recall that the inverse of Y represents the quantum lifetime of a particle in the 
corresponding resonant state escaping into the leads. Moreover we assume that the particles 
are uniformly distributed inside the sample and diffuse until they reach the boundaries, where 
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they are absorbed. Then we can associate the corresponding lifetimes with the time t R ~ 
1/r.R ~ R 2 /D a particle needs to reach the boundaries, when starting distance R away. This 
classical picture can be justified for all states with T > T c i ~ D/L 2 . The relative number of 
states that require a time t < t R in order to reach the boundaries (or equivalently the number 
of states with r > Tr) is 

v- m t(T R )= rr(T)dr~t>M (8) 

where S(t R ) is the area populated by all particles with lifetimes t < t R . In the case of open 
boundaries we get 

L 2 -(L- 2R) 2 [Y7i r c 



Pmt(r*) y — 2 '- ~ \ =-- — =— ■. (9) 

For > T c ; the first term in the above equation is the dominant one and thus Eq. JQ follows. 

Here it is interesting to point that a different way of opening the system might lead to a 
different power law behaviour for V(T). Such a situation can be realized if instead of opening 
the system at the boundaries we introduce "one-site" absorber (or one "lead") somewhere in 
the sample. In such a case we have 

V^ATn) ~ — - — - — (10) 

The above result is valid for any number M of "leads" such that the ratio M/L 2 scales as 
1/L 2 . In Fig. 2b we report 'Pint(r) for the case with nine "leads" attached somewhere to the 
2D sample. 

A straightforward generalization of our arguments for 3D systems in the metallic regime 
gives V int (T R ) ~ (T c i/T R f 5 -2(T cl /T R ) + (4/3) (r c/ /r fl ) 15 which for T R > T ct leads to 
the same universal expression as in Eq. (JJJ. Similarly, the analogue of Eq. in 3D is 

v int (r R ) ~ (r ci /r fl ) 15 . 

Wigner delay times distribution. Our theoretical understanding will be based on the fol- 
lowing relation 

L 2 

n— 1 v ' ' 

which connects the Wigner delay time and the poles of the S— matrix. Let us start with the 
far tails. It is evident that large times tw(w) ~ T^ 1 correspond to the cases when oj ~ oj n and 
r„ < 1. Then for the distribution of delay times we obtain V(r w ) ~ / dTV(r)TS(T W - 1/T). 
Using the small resonance width asymptotics given by Eq. JQ) we get the log-normal law of 
Eq. ©. 

Now we estimate the behaviour of V{tw) for tw < T^ 1 . In this regime many short-living 
resonances contribute to the sum Hll|) . We may therefore consider Tyy as a sum of many inde- 
pendent positive random variables each of the type T n = T n x n , where x n — Su>^ 2 . Assuming 
further that 5u> n are uniformly distributed random numbers we find that the distribution 
V(x n ) has the asymptotic power law behaviour 1/xf/ 2 - As a next step we find that the distri- 

3/2 3/2 

bution V{T n ) decays asymptotically as 1/t„ where we use that V(T n ) - 1/r,/ . Then the 
corresponding V(tw) is known to be a stable asymmetric Levy distribution L^\(tw) of index 
(j, = 1/2 [25] which has the form given in Eq. J3J) at the origin. We point out here that the 
asymptotic behaviour V(t\y) ~ ^/ T w 2 emerges also for chaotic/ ballistic systems where the 
assumption of uniformly distributed 5u n is the only crucial ingredient (see for example [2]). 
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Since Tw = 2i=i T q> we expect the behaviour of the distribution of proper delay times 
V(T q ) to be similar to V(tw) for large values of the arguments (for tw > 1 we have tw ~ 
r™ ax ). Moreover, from the numerical point of view V(T q ) can be studied in a better way 
because a larger set of data can be generated easily. Our numerical findings for V(T q ) are 
reported in Fig. 3 and are in nice agreement with Eq. (J2J), even for moderate values of r q . We 
stress here that the dashed lines in Fig. 3, have slopes equal to Cp taken from the corresponding 
log-normal tails of V(T). 

We thank R. Fleischmann and Y. Fyodorov for many useful discussions and comments and 
C. Beenakker for pointing a misprint in Eq. (0J. (T.K) thanks U. Smilansky for stimulating 
his interest in quantum scattering. This research was supported by a Grant from the GIF, 
the German-Israeli Foundation for Scientific Research and Development. 
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